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Abstract. Multi-valued logical models can be used to describe biologi¬ 
cal networks on a high level of abstraction based on the network structure 
and logical parameters capturing regulatory effects. Interestingly, the dy¬ 
namics of two distinct models need not necessarily be different, which 
might hint at either only non-functional characteristics distinguishing the 
models or at different possible implementations for the same behaviour. 
Here, we study the conditions allowing for such effects by analysing 
classes of dynamically equivalent models and both structurally maxi¬ 
mal and minimal representatives of such classes. Finally, we present an 
efficient algorithm that constructs a minimal representative of the re¬ 
spective class of a given multi-valued model. 
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1 Introduction 

Logical modelling has been used for describing and analysing biological systems 
such as gene regulatory networks for many years now, starting most notably with 
S. Kauffman’s work on Boolean networks The resulting models can be viewed 
as non-uniform cellular automaton [2] with varying diameter. The underlying 
idea of capturing a system in terms of network structure and logical rules that 
govern the change of components has later been extended by R. Thomas [TU] • In 
his multi-valued logical framework components are allowed to adopt more than 
two activity levels and so-called logical parameters determine the effect of sets 
of regulators on a target component. The update scheme yielding transitions in 
the state space of the model can be chosen to allow for either all indicated value 
changes to be changed in one step, i.e. synchronously, or for only one possible 
component update per step, i. e. asynchronously. In the asynchronous update, as 
introduced by R. Thomas, this results in a non-deterministic transition system 
that generally is harder to analyse than in the synchronous, deterministic, case, 
but usually includes more realistic trajectories la¬ 
in application, constructing a model, i.e., specifying the network structure 
and assigning the logical parameters, is a difficult task. Usually not enough 
information is available to make unambiguous choices. Consequently, comparison 
of transition systems of many models becomes necessary. In this context, it is 


an interesting observation that different models may give rise to the same state 
transition system. For Boolean networks this phenomenon has been investigated 
and is directly related to superfluous edges— edges without detectable dynamical 
effect—in the graph capturing the network structure (see e.g. [5]). In the multi¬ 
valued setting however, not only existence but also strength of a regulatory effect 
can be captured in the model. This allows for models differing both in structure 
and parameters to generate the same transition system, even if only functional 
edges are considered. This represent the full dynamical equivalence, as opposed 
to what is usually understood as dynamical equivalence in reduction studies, 
e.g. [7], where the focus is only on the stabilizing behaviour only. 

For a full understanding of this problem in the multi-valued case, the notion 
of functionality or observability needs to be extended from the straight-forward 
intuition as it is usually used (see e.g. [5] [3]). In this paper, we firstly present the 
proofs of results given originally in that allow to describe the equivalence class 
of models generating the same transition system using a structurally maximal 
representative. This already allows to pinpoint the key properties to be focused 
on, namely self-loops and strength of effects as modelled by the value of a given 
logical parameter. However, the approach is not suited to application due to 
the high complexity. Extending the results in [^, we then show in the second 
part of this paper that we can also systematically derive a structurally minimal 
representative of the equivalence class of a given model exploiting a resolved 
notion of observability and present an algorithm implementing our ideas. 

2 Background 

We start by introducing the relevant notions, with a simple illustrative example 
given in Fig. Hal 

A multi-valued regulatory graph (RG) is a triple G = {V, E, p) where: 

— F is a set of named components, 

— E C V X N+ X F is a set of regulations, 

— p : F —>■ N"*" assigns the maximal activity level to a component. 

It is required that for all (u, n,v) G E the non-zero integer n, called threshold, is 
in the range of the source component i.e. 0 < n < p{u). We also use an additional 
function 9 : V x V ^ 2”^^ which provides the thresholds of all edges between 
two vertices and is defined as 9{u,v) = {n | {u,n,v) G E} where u,v gV. Note 
that 9{u, ■(;) = 0 if there are no edges from u to v. 

We use the symbol Q to denote the set of multi-valued RGs. Also, in the 
following we use the notation Xi^-n for a substitution such that in a tuple x the 
value indexed by i is substituted for n. 

2.1 Discrete kinetic parameters 

We denote the state space of a system with S = Note that the 

state space is shared among the graphs that have the same function p and thus 
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Fig. 1: Toy example, a) A network G with the parametrization K. c) A reduced 
toy network G' with the parametrization K'. Note that for any x € [0,/9(u)] 
and y e [0,/;?(«)] it holds that F^{x,y) = F^'{x,y) and F^{x,y) = F^'{x,y). 
Therefore we have b) the transition system Ta{K) = Tq{K'). 


also the same V. We will use Op = {(V,E,p') | p = p'} to refer to the class of 
graphs that share the state space, e.g., the graphs in Figures ITal IT^ belong to 
the same class. 

The set S represents all the qualitatively different configurations of a system. 
However, each component is dependent only on the values of its regulators. An 
equivalence class on S w.r.t. regulation of a component u € H is called the 
regulatory context. To define the relevant notions we first describe the activity 
interval of a regulator. For formal reasons we consider an extended threshold 
function 0 with 0{u, v) = 0(u, v) U {0, p(u) + 1} for all u,v € V. Then 

G 0(u,v),j < k,-'(3l e 0(u,v)(j < I < fc))} (1) 

is the set of activity intervals of u in regulation of v. Here, the intuition is that 
the regulator effect of u on u is constant in each of the intervals of Note that 
[jly = [0, p(u)], even in the case that there are no edges from u to v. 

The set of regulatory contexts of v is then denoted and defined Gy = IIugv ■ 
For each v gV, a regulatory context tv G Gy is a product-of-sets where uju is the 
u-th set of the product for each u G V. For convenience we use comparisons of 
activity levels and intervals. We say that n < ujy, n G [0, p(u)] iff for each m G ujy 




























we have n < m. Also note that the intervals are not overlapping, therefore if we 
extend < to pair of intervals then < naturally forms a total ordering. 

The dynamics of the RG are given via integer values, called logical parameters, 
assigned to each context. We use a parametrization function Ky : fiy —>■ [0, p{v)] 
for each v G V. The parametrization of a regulatory network (V, E, p) G G is then 
a tuple K = {Ki, ..., K^y^). Lastly, we denote ICq the set of all parametrizations 
of the RG G, called the parametrization space. 

To denote whether an edge is superfluous in the system, we use the observabil¬ 
ity constraint. Intuitively we say that an edge is not observable if the parameter 
value is the same for all pairs of regulatory contexts that differ only in presence 
of the said edge. In Figure [T^ it is apparent that the change in the oJu does not 
cause a change in the parameter Ky(iu), i.e., the edge {u, I, v) is not observable. 

We now formalize this notion. In the rest of the article we will use for 
an edge {u,n,v) the symbol n_ G 0{u,v) for a value such that there is no 
m G 0{u,v) where n- < m < n, and G 0{u,v) for a value such that 
there is no m € 0{u,v) where n+ > m > n. Note that n_,n+ always exist as 
0, {p{u) + I) G 0{u,v) — 9{u,v). Now, we say that {u,n,v) is not observable in 
parametrization K G Kq iff for each uj G fiy such that = [n,n_|_) we have 

2.2 Asynchronous transition system 

Having a RG G = (V, E, p) and a parametrization K G ICg we can fully describe 
its dynamical behaviour as a transition system (TS) over its state space S. 
This is a directed graph (S', —>■) where —>■ C S x S is the transition relation. As 
mentioned, we are interested in asynchronous dynamics which means that the 
transition relation is non-deterministic. 

First K is converted into a so-called update function E^ = {Fff)y^v where 
Fff : S —>■ [0,p(n)] for all v G V. Here we exploit the fact that for each s G S 
and for each v G V there exists a context ui G fiy such that s G 
simplify the notation we will further write s G oj instead of s G YIugv^'^- 
every v G V we obtain the function F^ : S ^ S from a parametrization Ky as 

{ Sy 1, if Sy < Ky{uj), S G UJ, 

Sy, \i Sy = Ky{uj),S G IZJ, (2) 

Sy — I, if Sy > Ky(uj), S G UJ. 

Note that we have only three options how the value of a component can change, 
namely either increase by one, remain constant, or decrease by one. This provides 
us with a certain notion of direction of a derivative (positive, zero, or negative). 
In the following we will be using the term S{Fy){s) = Fy{s) — Sy to describe the 
partial derivative in the dimension n in s G S'. 

Having F^ , we now assign each parametrized RG a TS via the function 
Tg '. ICg {(S, —>■)} where Tg{K) = (S, —)•) such that: 

Vs, s' G S : s ^ s' (s ^ s') A {F^ (s) = n) A (3n gV : s' = Sy^n)- (3) 
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Fig. 2: Example of uncertain reduction. The network b) is in a non-canonical 
form. The canonization a) however makes the edge (r;,l,w) observable, even 
though it is superfluous as illustrated by c). 


In the following, we will compare the resulting TSs generated by RGs with 
the same state space, ie., those in Qp for some p. We denote this set of TSs 
Tp = {(5', —>•) I Tc{K) = (S', —>■), G € Qp, K € /Cg}- An example of a shared TS 
is given in Figure [TJ 

3 Conservative graph manipulations 

We now investigate the cases where different parametrizations of a single graph 
generate the same TS. From (jS]) it is clear that two functions ^ will lead 
to distinct TSs, while coinciding functions F^, F^ lead to the same dynamics. 
We therefore focus on describing the situations where for K ^ K' we still have 
F^ = F^'. 

Consider the simple example in Fig.[2j We see that putting 1}, {1, 2}) = 

1 instead K{{{), 1}, {1, 2}) = 2 still yields F^ (0,1) = F^{Q, 1) = 1 and there¬ 
fore the TS remains the same. This illustrates that, other than in the Boolean 
case, information on parameter values may get lost when deriving the update 
function. In Fig. [23 we have a case where the parameter value lies outside its 
context and by incremental change we leave the context even before the value 
can be attained. 

We now define the notion of a canonical parametrization that prohibits such 
effects. Observe that a value change in v can cause the change of context only if 























V regulates itself. Therefore we say that K € K-g is canonical if and only if 


Vu e V, Vw e f2v,cav = [j, k) : {Ky{uj) > j — 1) A {Ky{oj) < k). (4) 

We also denote Cq Q ICq the subset of canonical parametrizations in JCg- 

We can obtain a clear correspondence between K and if all the contexts 
contain just a single state, so that no ambiguities are introduced in ([2]). This 
partition is achieved when only considering complete graphs. For clarity we add 
that {V, E,p) G Q is complete if and only if for all u,v gV and every n G [1, p{u)] 
the edge (u, n, v) is in E. This gives us the following theorem: 

Theorem 1. For each G = {V,E,p) G Q it holds that if G is complete then Tq 
defines a bijection between Cg and Tp. 

Proof. Let Tg{K) = (5, —>■) for some complete G = {V, E,p) G G such that K G 
Cg- The one-to-one correspondence between F^ and {S, —>■) immediately follows 
from We therefore need to show that there is also such a correspondence 
between K and F^. 

First, it is important to note that if the graph is complete, each regulatory 
context depicts only a single configuration. This is because if G is complete then 
by (HD we have: 

\lvGV:Gy=W_ 1)’ [1> 2),..., [p{u),p{u) + 1)}. 

u^V 

Since each context is a singleton, each component has its value fixed. In that 
case, canonicity requires that the parameter in the context differs from the value 
only by 1. Therefore we have only three options for the parameter value. More 
precisely, by substituting we have 

G V,ys G S : St, — 1 < Kv{{s}) < s„ -|- 1. 

Then in such a case, can be written as 

{ Sy 1; if Sy “hi - 

s„, if S.„ = iF„({s}), 

Sy 1, If 1 — 

from which we immediately get Ftf-(s) = Ky{{s}). Thus ([2D corresponds to: 
yv GV,\/s G S : (iF„({s}) = n) A (iF„({s}) ^ s„) s s„<_n. 


□ 

Based on this theorem, we can consider a complete graph with canonical para- 
metrization as a representative of a class of models with the same behaviour. 
Now we show that it is possible to convert any graph with some parametriza- 
tion into a complete graph with canonical parametrization, while keeping the 
dynamics unchanged. 


First, we focus on the canonization function Can : {((?, iC) | G G G,K G 
ICc} —>■ {(G, C) I G G Q,C G Cg}- For each component v G V and for each 
regulatory context to G Gy with ujy = [j, k) we construct G as follows: 

fj-l, if < j - 1 

Gy{uj) = ■!/(:, if Ky{uj) > k 

\^Ky{uj), otherwise 

The goal is to avoid that the parameter value cannot be reached in one 
transition from any state in the context. The procedure is illustrated in the 
conversion from the parametrization K in Figure I2bl to the parametrization 
K' in Figure [2al We now prove that this procedure indeed yields a canonical 
parametrization for any RG that shares the TS with the original one. 

Lemma 1. Canonization is correct. For all G G G and all K G /Cg it holds that 
if Can{G, K) = (G, G) then G is canonical. 

Proof. There are two options for K not to be canonical. The first option is that 
3v GV,3 uj G Gy, ujy = [j, k) : Ky{uj) <j-l 
but then G^(a;) = j — 1, so G is canonical. The second case 
G V,3u! G Gy,ujy = [j, k) : Kyito) > k 

can be treated analogously. □ 

Lemma 2. Canonization is conservative. For all G G G and all K G /Cg it 
holds that if Can{G, K) = (G, G) then Tg{K) = Tg(G). 

Proof. Recall that the TSs Tg{K) and Tg{G) are fully defined by F^ and F'^, 
respectively. We therefore need to show that F^ = F'^. 

For all u G R and for all oj G Gy the value Gy{uj) is set based on one of the 
three following cases. 

First consider the case that Ky{uj) < ojy. Denote ujy = [j, k). For all s G w 
it holds that Sy > j — 1 and then Cy{(jj) = j — 1. This means that for all s G w 
both Ky(uj) and Gy {tv) are smaller than Sy and therefore for each s G w we have 
= sy-l = F^{s). 

The case that Ky{uiy) > ujy can be treated analogously. 

The third case is that we have Ky{ujy) = Cy(ujy) and thus by definition 
F^(s) = for any s G . □ 

Corollary 1. For any K holds: If j = Ky(ijj) < ujy (resp. j = Kyito) > tOy) 
then replacing j with f, 0 < j' < tOy (resp. tOy < j' < piv)) is conservative. 

Second, we extend the structure of a graph using the completion function Comp : 
{(G, iC) \ G G G, K G /Cg} ^ {{G,K) \ G G G, K G /Cq}. If G is complete, we 
map {G,K) to itself. For an incomplete G = {V,E,p) and some K G /Cg we 


consider the non-empty set of missing edges E = {(u,n,v) \ u,v S V^n G 
[l,p('«)], {u,n,v) ^ E}. Assume that the set of all possible edges has some or¬ 
dering. We extend the graph G to G such that G = {V,E\J {min{E)}, p). As an 
example you can see that the RGs in Figure [T^ [T^ differ from each other only 
by the presence of a non-observable edge. Likewise for the GRs in Figure [?bl[2cl 
To extend the parametrization K to the new structure we first observe that G 
gives rise to new contexts that were obtained by partitioning some context of G 
into two. To preserve the dynamical behaviour we simply assign the parameter 
value of the original context to both resulting new contexts. Have (u, n, v) = 
min{E). For each v €V and for each ui € f2y we then create K as: 

- J A:«(w) a V V W (uJu = [j, k) n+)), 

otherwise. 

We now prove that for an incomplete RG we can use the completion procedure 
to add a new edge while retaining the dynamics. 

Lemma 3. Completion is sound. For all G € G and all K G ICg it holds that if 
Comp{G, K) = (G, K) then G € G and K G Kq. 

Proof. We have u,v € V and n G p{u), therefore by definition of G we have that 
if {V, E,p) G G then (V, E U {(u, n, u)}, p) G G- 

The variables n_,n+ exist since 0 is always a possible choice for n_ and 
p{u) + 1 for n+. For any m G 0{u, v) we know that (0 < m < p{u) + 1). 

From © we know that the only change occurs in the interval We therefore 
only need to show that Ky is extended to the affected contexts. Since n_,n+ 
exist, for any uj G f2y we have that uju S {G, • ■ • > [n-jn+), ..., ik}- Therefore for 
any w G we also have S Therefore Ky is defined on the whole 

f2v for each v gV. □ 

Lemma 4. Completion is conservative. For each G G G and for each K G Kg 
it holds that if Comp{G, K) = {G,K) then Tg{K) = Tq{K). 

Proof. We have that K differs from K only in a context Cj G Qy with Cj = [j, k) 
where either j = n_ or fc = n+. 

Assume there is some s G w for which F^{s) ^ F^{s). But we know that 
[j, k) C [n_,n+) and therefore s G This would however imply that 

also Kyifby) ^ Ky(u)u^_^^n-,n+}), which contradicts the definition of K. We there¬ 
fore have that F^(s) = F^{s). □ 

Since the completion procedure adds only one edge at a time, we need to repeat 
the procedure. This is captured in the following Lemma. 

Lemma 5. For G G G, K G ICg, consider the recursive sequence Comp(G, K), 
Comp{Comp{G, K)),... 

This sequence converges to a fixed point (G°, K'^) and G° is complete. 


Proof. The set E of missing edges in G is finite as V is finite. For each v € V 
also [1, /3(f)] is finite. In each iterative use of Comp the size of E is decremented 
by one. The recursive sequence becomes constant when E is empty, indicating a 
fixed point of Comp. By definition, is complete. □ 

Combining all the statements above, we arrive at our final theorem; 

Theorem 2. Let G,G' &Q, K € Kg, K' G Kq’ and denote Comp*{G, K), resp. 
Comp* (G', K') the fixed points derived from iterating Comp starting in {G,K) 
resp. {G',K'). 

Then Tg{K) = r^(iC') iff Can{Comp*{G, K)) = Can{Comp*{G',K')). 

Proof. We now know that Can(Comp*(G', K)) and Can(Comp*(G', K')) are canon¬ 
ical and complete. The equivalence follows from Tq being a bijection, as proven 
in Theorem [TJ □ 

4 Network minimization 

Canonization and completion provide mathematical insights into the problem 
of dynamical equivalence. However in application it is more useful to have a 
minimal, rather than maximal structure. Also, while the canonization provides 
a good intuition about what the actual behaviour for each context is, it can have 
side-effects like converting a non-observable edge to observable, as is illustrated 
in Figure [D We therefore introduce another form of parametrization, named 
normalized parametrization, which prevents such effects but is more involved. 
Using a normalization procedure we then obtain a parametrization which is 
amenable to minimization. This section is divided into five consecutive steps: 

1. We introduce a notion of observability in the TS which allows us to see 
whether an edge is observable based on the transitions in the TS. 

2. We introduce a notion of a monotone target value (MTV) of a component. 
This value keeps the observability properties of TS, but is shared for a whole 
context. 

3. We show how to compute the MTV from a parametrization and consequently 
how to compute a normalized parametrization. 

4. We show that every edge that is not observable in the TS is not observable 
in the respective normalized parametrization. 

5. We introduce the minimization function for RGs based on normalized para¬ 
metrization and explain how to test equivalence via minimization. 

4.1 Observability in transition systems 

We have already defined the notion of observability in the parametrization. How¬ 
ever we are mostly interested in the observability since it has implications on 
the dynamics. We now show how the it can be evaluated in the TS. 

Intuitively, for an edge to be not observable there must be only a single value 
towards which the component evolves, no matter whether the regulator is above 


the thresholds of the said edge or below it. Formally for an RG G = (F, E, p) is 
an edge (u, n,v) £ E not observable in {S, -^) = Tc{K) iff: 

Vs € S,Su € [n_,n+),3fc G [0,p(w)],Vj G [n_,n+) : 

6{Ff^){su^j) = Sgn{k - {su^j)v) , (5) 

where Sgn : Z —>• {+1, 0, —1} is the usual sign function. 

Consider the example in Figure[TJ It is easy to see that the regulation (m, 1, v) 
does not have any effect, since the left and right half of the TS are identical. 
Take in particular the example of the state {v,u) = (0,0). We choose k = 2. 
Then for j = 1 it holds that S{E^){0, !) = +! = Sgn(2 — 0) and for j = 0 it 
holds that (5(F)f )(0,0) = +1 = Sgn(2 — 0). 

From the definition of the function F and its derivative we easily see that: 

Ws € S,Sy eoj £ Hy : S{F^){s) = Sgn{Ky{u}) - Sy). (6) 

This illustrates that the TS non-observability is the kind that we are interested 
in, since it relies on the actual dynamics of the network as captured in F^. It 
is a stronger notion than the corresponding parametrization based one, since 
observability in the TS implies observability in the parametrization. We show 
the contraposition of this statement in the following lemma: 

Lemma 6. Have {u,n,v) £ E not observable in K £ K-g- Then (u,n,v) is not 
observable in Tg{K). 

Proof. By definition of observability (Section 12.11) we have that for each oj £ f2y 
such that ujy = [n,n+) it holds that Ky(uj) = R'«(wu<_[„__„)). We therefore can 
set k = Kyiio) = and from ([6]) we immediately see that ([5]) is 

satisfied for the whole range [n_,n+). □ 

4.2 Monotone target value 

To relate the dynamics captured in a TS with a parametrization value, we in¬ 
troduce the notion of monotone target value. Intuitively, for a state s G S' an 
MTV is a value towards which the component value Sy evolves if we traverse 
only in the dimension of v until either stabilizes or an opposite effect takes 
place. This idea is strongly linked to the derivative of the update function. 

Consider the example in Figure fTbl Under the influence of the edge (u, 1, v) we 
see the trace (0,1) —>■ (1,1) ££ (2,1). Here the MTV for component v in the state 
iv,u) = (0,0) is 2 since <5(Ff)(0,l) = +I = 5(Fj^(l,I)) ^ <5(Ff (2,1)) = -I, 
he., at the level St, = 2 an opposing effect takes place. 

For a Tg(K) = (S, —>■), v £ V, u) £ Gy and any s G w we denote the MTV 
by (Ff )’"°”(s) defined as: 

isy, if (5(F„^)(s) = 0, 

{F^r°^{s) = I rmn{j > s„|^(Ff )(s„^,) ^ +1}, if <5(F„^)(s) = +1, (7) 
[max{j < Sy\S{F^){sy^j) -I}, if 6{F^){s) = -I. 


Also note that: 


'iseSyveV: S{Fi^)is) = Sgn{iFi^r°^{s) - s,). ( 8 ) 

The MTV can now be related to the observability in the TS. In particular, we 
can rewrite m as: 

Vs € S,Su € [n_,n+),Vj G [n_,n+) : 

5{F^){s^^,) = Sgn((F„^)™°"(s) - (s„^,),) (9) 


Lemma 7. The non-observability conditions ^ and m are equivalent. 

Proof. Clearly, ([H]) implies ([5]) since we can set k = (T'J^)™°"’(s). For the other 
direction we show that if ([9]) does not hold, then ([5]) cannot hold either. 

If ([9]) does not hold, then there is a state s G S' and some j G [n_,n+) such 
that S(FJ^)(s') 7 ^ Sgn((Fj^)"*°"(s) — s^) where s' = Su<_j. We distinguish three 
cases based on the value of S(Fj^)(s') and then again three cases based on the 
difference between Sy and sj,: 

Case S(Ff^)(s') = 0 and Sgn((Fj^)'"°”(s) — s^) ^ 0: 

— If s„ = Sy, then for any k certainly Sgn(/c — s^) = Sgn(fc — Sy). Also 
since Sgn((Fjf^)’”°"’(s) — sj,) 7 ^ 0 we get (Fi^)™°"(s) ^ Sy = Sy, and thus 
(5(Fjf(s)) 7 ^ 0 according to ([ 8 ]). Then if there is a fc s.t. S{Ff^){s') = 
Sgn(fc — Sy) then also Sgn(fc — Sy) 7 ^ 6{F^{s)) and therefore (O does not 
hold. 

— If s^ > Sy, it follows that u = v. From (O we get {Ff^)'^°"‘{.s) > Sy since 

(5(Fjf )(s') = 0 indicating an effect change in the only problematic case that 
6{Ff^)ls) = —1. Since Sgn{{Ff^)"^°^{s) — s'y) ^ 0 we have strict inequality 
((Fjf-)"''°"'(s)) > s^. Subsequently by ([7]) there exists I such that Sy > I > Sy 
and 6{F^){sy^-i) > 0. We have that Sy G [n_,n+) and S{F^){s') = 0, so 
to fulfil (O the k must be chosen as k = s'y. But 5{F:lf){sy^i) > 0 and 
I G (s^,s«] C [n_,n+) and therefore k needs to satisfy k > = 1. 

Together we get k = s'y < I < k which is a contradiction. 

— The case Sy < s'y can be treated analogously to Sy > s'y. 

Case 6{Ff^){s') = +1 and Sgn((Fjf')"*°"(s) — s'y) < 0: 

— If = s'y, then, since Sgn((F,f-)’”°"’(s) — s^) < 0 and thus < 

s'y = Sy, we have S{Ff^){s) < 0. Then /c in ([SJ needs to satisfy k < Sy. Also 
S{Ff^){s') = +1 = Sgn(fc — s^) and thus k needs to satisfy k > s'y leading to 
a contradiction. 

— The case Sy > s'y is impossible, since Sy> s'y > (F^)’"°”(s), so 5{F^){s) = 
— 1 and therefore by © also 5{F^){s') = — 1 which is a contradiction. 

— If S.U < s'y, then by © there is I such that Sy < I < s'y and 6{Ff^){sy.f-i) < 0 . 
As I G [n_,n+), a k satisfying condition ® must be such that Sgn(fc — 
{sy^-i)y) < 0 and thus k < 1. But S{Ff^){s') = +1 and a suitable k must also 
satisfy s'y < k. Together we again have the contradiction k < I < s'y < k. 


Case S{F^){s') = -1 and Sgn((i^^)'"°"(s) - s^) > 0: 
This case can be treated like the previous one. 


□ 


4.3 Normalization algorithm 

We can see that the MTVs allow for a straightforward test of ob¬ 

servability in the TS. Obtaining {F^)'^°'^{s) is however quite tedious. The size 
of S is exponential w.r.t. the set V and we have to unfold the TS to find the 
monotone paths characterizing the MTVs. In this section we show that all states 
of a context share their MTV and additionally that we can obtain the MTV from 
a context directly. 

We introduce the normalization function, described in Algorithm [1] that 
computes for each component and for each regulatory context of that component 
the MTV shared between the states of the context. 


Algorithm 1 Calculate Norm(A, n,a;) where 0{v,v) = {no, ...,nfc}. 

1 . = uJv 

2: if Ktj{u!) £ ljJv then 

3: Norm (A, a, a;) = (a;) 

4: else if A„(tj) < then 

6: if Kv{uj') > Ui — 1 then 

7: Norm(A, a,aj) = Ui — 1 

8: else 

9: Norm(A, a, a;) = Norm(A, u, ca') 

10: end if 

11: else 

12 . LO 

13: if Kv{ui') < rii+i then 

14: Norm(A, v,a;) = rii+i 

15: else 

16: Norm(A, a, w) = Norm(A, w, ca') 

17: end if 

18: end if 


In the algorithm we traverse through the contexts, rather than through 
states of a system, when looking for a monotone trace. As an example con¬ 
sider the A{ in Figure [2a1 and ca = ([0,1), [1, 2)). Then Norm{K',v,uj) = 
iVorm(A', w,= 2. Note that this coincides 
with the value in Figure [2H which actually has a normalized parametrization. 

The correctness of the approach is quite intuitive since a regulatory context 
is a subspace of the state space with uniquely determined target value. This 
means that either the behaviour is monotone or there is exactly one stable state 
which breaks monotonicity in both directions. Considering the definitions of 
(^A')mon(g) |.]gg derivative. Easy calculations for the three cases Ky{uj) < ujy, 






Ky{uj) > ujy and Ky(uj) € u>y for a context w immediately prove the following 
lemma. 


Lemma 8. For anyv G V and any u € fly we have 
for all s, s' G w. 

Due to this Lemma we can extend the notion of MTV to regulatory contexts 
u! SO that (Fjf for any s G w. Having this extension, we 

now prove the correctness of Algorithm [TJ 

Theorem 3. For v gV, uj € f2y it holds that (Ajf )™°"(a;) = Norm{K,v,u!). 

Proof. Linking back to the importance of self-regulation already seen in the 
canonization, we lead the proof by induction w.r.t. the distance in number of 
activity intervals of self-regulation of v between the context and its MTV. We 
want to show that Norm(iG, v, uj) returns the correct value after at most as many 
recursive calls as is the distance between the context and its MTV. The notion 
of distance needed for this is defined as a function —>■ N with 

(^) = Max(|{A G If \ lOy < A < A^}|, |{A G If \ ujy > A> 

where Ayj G If is the activity interval where it holds that G A^. 

Now we prove the theorem by the means of induction. Note that the proof 
does not constitute an invariant of the algorithm, as it proceeds in the other 
direction than the algorithm itself. In particular, we start by showing that for 
all the contexts that have their MTV within them or on their boundaries, the 
algorithm ends immediately with the correct value. Then we proceed to show 
that if the normalized parameter of any context whose distance is m is correct 
and known, then the normalized parameter of a context whose distance is m -I-1 
can be correctly determined by calling Algorithm [T] once. 

Base of induction (distance 0 and 1): 

If Dy (lo) = 0 then we know that G ojy. This implies that there is 

a state s Goj such that Ky{w) = Sy and S{Ff^){s) = 0 according to ([2]) and ([T]). 
Therefore = Ky(uj) = Norm(Ar, n,w), as set on the lines 2, 3. The 

recursion depth is 0. 

If Dy{uj) = 1, denote w' = (jJy.^A such that G A. It follows from 

Dy{uj) = 1 that Dy{uj') < 1. We distinguish the two options: 

— Dff{uj') = 0: Then the above argument repeats and there is s G w such 

that Kyiuj') = Sy = (Fif-)"‘°"(w'). By the definition of the MTV we get 
(Fjf-)™°”(w') = (F^)™°”(a;). Based on the ordering of uj,uj' we arrive either 
on the line 9 or 16 of the algorithm and state correctly that = 

Norm(Ar, = Norm(Ar, u, w') = s„. The recursion depth is 1. 

— Df'(uj') = 1: Since Dy{uj) = 1 it follows that (F^)'"°"(a;') G ujy. In this 
case the respective MTVs take the adjacent values of the boundary between 
LOy and cu'y. In the case that uj > uj' we arrive on line 7 in the algorithm 
and assign Norm(Ar, z;, w) = n — 1. This is correct as in any state of w we 


monotonously update towards lu' and as we enter uj' by crossing the boundary 
value n, we change the direction back towards w, breaking the monotonicity. 
Analogously the correct value is assigned for the case uj < uj'. The recursion 
depth is 0. 

Induction step (distance over 1): 

The induction assumption is that in at most recursion depth m > 1 the value of 
any uj' G such that {oj') < m is correctly set and consider now (uj) = 
m + 1. 

In case uJv > since {uj) > 1, we have an uj' such that 

D^{uj') = TO > 1 and uJv > uj'^ > From the definition of the MTV it 

follows that all states in both uj and uj' monotonously decrease under , there¬ 
fore also uj', > (Ff which gives us (Ff = (Ff In the 

algorithm this is assured on line 9 and by induction hypothesis (F^)™°"(w') is 
correctly determined by the algorithm in at most to recursions, giving us the de¬ 
sired result for (F^)’"°”(a;). The case that uj, < (F^f-is again analogous, 
leading to line 14 instead of 9. The recursion depth is now to -I- 1. 

Since the set of intervals is finite and the recursion traverses monotonously, 
we terminate in the recursion depth of at most Max({p(u)|u G V})- □ 

Using the normalization function we can, similarly to canonization, create a 
conservative and sound transformer on parametrizations. We extend Norm to 
a function Norm : {(G,F) | G G C?, F G K-g) ^ {(G,iV) | G G e,iV G Mg) 
where Mg Q f^G is the set of normalized parametrizations of G = (U, F, p) and 
Norm(G, K) = (G, N) where N is defined by 

Vu G U, Vw G Gj, : N,{uj) = Norm(F, u, w). 

We have proven correctness of normalization already in Theorem [31 so it only 
remains to prove that normalization is conservative. 

Lemma 9. Normalization is conservative. For all G G Q and every K G Kg it 
holds that if Norm{G, K) = (G,N) then Tg{K) = Tg{N). 

Proof. Observe that if K,{uj) G ujy then Norm(F, w,a;) = K,(uj). In Corollary[T] 
we have shown that if K,(uj) < uiy, then it is conservative to replace K,(uj) with 
any I G N such that I < uJv, which is also the case in Algorithm |TJ The same 
holds for the case that Ky{uj) > ujy. □ 

4.4 Observability in normalized parametrization 

We have seen now that observability in the sense of an actual dynamical effect 
should not be evaluated based on the parametrization but rather on the TS. All 
information needed to construct a TS is captured in the MTVs due to its relation 
to the derivative and thus the update function. At the same time, an important 
aspect of parametrizations is shared, namely that the MTV stays fixed within 
a context. This allows us to link observability in parametrization and TS, as is 
shown in the following theorem that complements Lemma 6. 


Theorem 4. For all G G Q and every K G K-g it holds that if Norm{G, K) = 
{G,N) then every edge that is not observable in Tg{K) is not observable in N. 

Proof. Assume that the above does not hold, i.e. there exists an edge {u, n, v) G 
E s.t. dH) holds, but also it holds that: 

3cU G , UJu — [^5 — ^ui—[n_ ,n) ■ (^) (^'^) ■ 

Case u ^ v: 

Note that in this case we have a;„ = 

We have S{Fff){s) = d{Ff^){s) for all states s G uj as can be easily deduced 
from Lemma|3 For all s e w, s' € w*- we have (F^)"''°"'(s) = Ny{uj) ^ Ny{uj^) = 
(F'J^)”''°"(s'). It follows that we can only meet the condition Sgn((Fjf’)”''°"(s) — 
Sy) = Sgn((F„^)”"°"-(s) - s^) = 6{F^){s') for all s e a;,s' = Sy^jJ G w'*' as 
demanded in if and only if either Ny{uj) < cuy A Ny{uj^) < or Ny{uj) > 

LOy A Ny{uj^) > UJ^. 

First consider that {Ny{uj) < ojy) A {Ny{ijj^) < ujf). If the condition on the 
line 7 is satisfied for both oj and w'*' we immediately see that Ny{uj) = Ny{uji'). 
If it is satisfied for exactly one, then apparently we break the requirement ® 
as for s Gujy we have (s„^„_)^ = s„ but 5{F^){su^n.) ^ 5{F^){s). 

We therefore meet the condition on the line 8 and from the line 9 we know that 
for uj' as defined there Ny{uj') = Ny(uj) ^ Ny{uj^) = Ny{{ijj^)'). Since Ny{uj') ^ 
Ny{{uj^y) we have again that Ny{oj') ^ uj'y and the same for (w*-)'. Therefore 
it again must hold that Ny{uj') < ujy and A^„((wl')') < Apparently, the 

argument is recursive, requiring that for each ui' G f2y such that iVy < ujy it 
holds that Ny{uj') < ojy. But then ultimately Ny{ijj) < 0 which contradicts the 
definition of K. 

For the case that Ny^tv) > ujy a similar argument holds using the upper 
boundary Nyiio) < p{v). 

Case u = v: 

First note that in this case we have w' = for uj' as defined in the algorithm. 
In case that k = < n for any s G oj we execute Norm(Ar, u,w) in the 

algorithm and the condition on the line 4 is met. Then depending on the value 
Kiev'): 

— K{uj') > n — 1: then k = n — 1, meaning that for any s' G uj' we have 
(Fjf-= n — 1, otherwise the edge would be observable. But then 
also Ny{uj) = Ny{uj') which contradicts (US. 

— Kiuj') < n — 1: then according to the line 9 we have Norm(Ar, w,a;) = 
Norm(Ar, ?;,a;'), again contradicting (fTUl) . 

The case that k > n can be treated similarly. □ 

4.5 Minimizing the model 

Since we know that after normalization, all non-observable edges can be directly 
detected, constructing a reduction algorithm is rather straight-forward. Our min¬ 
imization process is closely related to the completion process in Section [3l 


We use the minimization function Min : {{G,N) \ G £ G,N £ Mg} —t 
{(G, iV) \ G £ G, N £ Mq} to eliminate the non-observable edges. If iG,N) is 
minimal, i.e. there are no non-observable edges, we map it to itself. Otherwise 
have E the set of non-observable edges in G = (V, E, p) and an arbitrary total 
ordering on E and (u,n,v) = min{E), then: 

1. G = {V,E\ {{u,n,v)},p), 

2. ,n_j_) ) [n_ ,n) ) ^v{^u-i^[n,nj^)'} • 

The nature of MTVs ensures that N is again in Mq. As an example consider the 
network in Figure [2bl The edge (u, Iju) is apparently not observable. We there¬ 
fore remove it from E and with that we set fVt,([0, 3), [0,1)) = A^u([0,1), [0,1)) = 
iV,([l, 3), [0,1)) = 0 and iV,([0, 3), [1, 2)) = iV,([0,1), [1, 2)) = iV.aij 3), [1, 2)) = 
2. Note that this coincides with the network in Figurewhich is in fact mini¬ 
mized and normalized. 

This procedure can be seen as an inversion to the completion as defined in 
Section[31 where we created two new contexts by splitting one, keeping the values, 
whereas here we merge two contexts with the same value into one. Note that 
the fixed points of Comp and Min are not dependent on the order on 1/ x N x F. 
We can therefore take arbitrary, but fixed, order and execute both Comp and 
Min according to this order. Then for Min(G,Af) = {G,N) with N £ M and 
G ^ G we have Comp(G,iV) = {G,N). This can be easily verified by applying 
the two operations successively. Since completion is sound and conservative, 
minimization in such corresponding cases is also sound and conservative. The 
only remaining case is that Min(G, N) = (G, N), but this is obviously sound and 
conservative too. 

Iteration of the Min function will then lead to a minimal structure, as demon¬ 
strated by the following lemma: 

Lemma 10. For G £ G, K £ ICq and Norm{G, K) = {G,N) consider the 
recursive sequence Min{G, N), Min{Min{G, N)),... 

This sequenee converges to a fixed point (G™, A^"*) and G™ is minimal, i.e. there 
are no other {G',K') s.t. Tq{K) = Tg>{K') and G' has less edges than G™. 

Proof. Since we change N only if we remove an edge and the edge set is finite, 
the existence of a fixed point is trivial. If G™ was not minimal, then inevitably 
there would have to be an edge which is not observable in Tq{K) but observable 
in N which contradicts Theorem |4l □ 

Now we can conclude the section with a theorem about equivalence checking 
through minimization, complementing the result of Section [3] and showing that 
the set of dynamically equivalent RGs has both maximal and minimal elements 
w.r.t. set inclusion on the regulators. 

Theorem 5. Let G,G' £ G, K £ K-g, K' £ KLg' and denote Min*{Norm{G, K)) 
and Min* {Norm{G', K')) the fixed points derived from iterating Min starting in 
Norm{G,K) resp. Norm{G',K'). 

Then Tg{K) = T^{K') iff Min*{Norm(G, K)) = Min*{Norm(G' ,K')). 


Proof. As in the proof of Theorem[5]we now know that both Min*(Norm(G', K)) 
and Min*(Norm(G', if')) are normalized and minimal. Due to Theorem 2] we 
know that in the minimization exactly the non-observable edges in the TS are 
removed. So if the two TSs are equivalent, the set of remaining edges is the same 
in both G and G'. 

Additionally we know that the parameter values are not changed during the 
minimization, only in the normalization where the value is set to the MTV of 
the states covered by the context. The only way how we could have a difference 
between Min*(Norm(G, AT)) and Min*(Norm(G', AT')) is that there is a state s G 
Tc{K) and a component v such that (A^)™°"(s) ^ But then by 

([7]) we have some r G Tc{K) such that 5{F^){r) ^ 5{F^'){r) and To{K) ^ 
Tg'{.K'). □ 


5 Complexity 

It has been shown that testing equivalence for two expressions over n variables 
is co-NP complete [T]. Since in a Boolean network a component v G V can 
be regulated by up to \V\ nodes, the update function can be an expression 
over \V\ variables and deciding whether two parametrizations are equivalent is 
therefore necessarily exponential w.r.t. the number of components. However in 
the minimization approach we can parametrize by the in-degree. For a given RG 
G = (V,E,p), denote v~ = {{u,n,v)\3n G N,3u G V,{u,n,v) G F} and set 
k = Max({|u“| I V G V}), then: 

Theorem 6. TIME{Min*{Norm{G,K))) = 0{\V\ ■ 2'= • k). 

Proof. There are at most 2^ different regulatory contexts for any u S V. As 
we call Norm(Ar, u,w) for each oj G Gy no recursion is needed—either we set 
the value directly or we set it equal to some other value that will eventually be 
known. We therefore call Norm(Ar, u, w) at most 2^ times for each v gV. 

In the minimization part we need to obtain the set E. For an edge (u, n, v) G 
E to detect whether {u^n,v) G E we need to consider all cu G Gy with = 
[n, n+) and compare them to the respective At most we need to do 

(|l7.i,|/2) G 0(2^) pair-wise comparisons. Since we have to consider at most k 
regulators of \V\ components, we obtain in total G(|R| • 2* • k). To remove the 
edge we remove the set of tested contexts, which again leads to G(| V| ■ 2^ • fc). □ 

Intuitively, since we only change the values in place or remove them, no additional 
space is needed. The spatial complexity is thus equal to the size of the input. 

6 Changing the npdate scheme 

Throughout the article we have been considering only the asynchronous update 
function, which was advantageous for the definition of non-observability in the 
TS and the monotone update. However we can extend our result to arbitrary 


update schemes where Fy{s) can be applied for any v €V and s G S, see [3] for 
examples. As an example, we consider the synchronous update. Then Tg{K) = 
(5, such that: 

Vs, s' G S : s ^ s' \/v GV : Fy{s) = s'y 

We can therefore immediately see that two TSs are again identical iff the respec¬ 
tive update functions are identical. Since the conversion from parametrization 
to update function remains unchanged, we have that two parametrizations are 
dynamically equivalent under the asynchronous update scheme iff they are dy¬ 
namically equivalent under any other scheme that depends on F. 

7 Conclusion 

In the article we have thoroughly investigated the notion of observability for 
multi-valued logical models of regulatory networks with an arbitrary update 
function. We have shown that, unlike in the case of Boolean networks, a mere 
test on equality of parameters in contexts with and without an edge is not suffi¬ 
cient to determine if the edge has a detectable impact on the network dynamics. 
As illustrated, at the heart of such an effect lies the connection of self-regulation 
and strength of regulation in general, as encoded in the parameter values. Here, 
self-regulations can act as amplifiers for the effect of other regulators. Conse¬ 
quently the same effect could be achieved without such a self-regulation simply 
by strengthening the original regulatory effect, i.e. by adapting the correspond¬ 
ing parameter value. 

One way to obtain an unambiguous representation of a given model that 
characterizes the whole class of dynamically equivalent models utilizes the no¬ 
tion of canonical parametrization. Here, in essence, the effects of regulations are 
standardized, but to do so, the underlying structure, i.e., the set of regulators 
for each component, has to be blown up. Although fully answering the mathe¬ 
matical problem of characterizing the equivalence classes, this approach suffers 
from high computational complexity in practice. 

Exploiting a refined notion of edge observability directly tied to the dynam¬ 
ical behaviour as captured in the transition system, we were able to construct 
a minimal parametrized RG w.r.t. the set inclusion on the set of regulators as 
representative of a dynamical equivalence class. The key step to make the con¬ 
struction feasible for application was finding a parametrization where the refined 
and the original notion of observability coincides. We presented an efficient and 
simple algorithm that allows us to obtain this normalized parametrization from 
an arbitrary one. 

Our thorough mathematical investigation uncovered a minimal representa¬ 
tive that is helpful not only in model comparison when modelling under uncer¬ 
tainties but also for understanding the regulatory effects underlying an observed 
behaviour. However, it should be mentioned that other models in the class might 
also carry interesting information. In essence they represent the different possi¬ 
ble mechanisms that can be used to implement a desired effect. In particular. 


in biological applications comparison between the actual and the possible im¬ 
plementations might help to gain a deeper understanding about the constraints 
governing biological systems. 

From the practical perspective, since the algorithm has very low complexity, 
we would suggest employing it in tools that focus on parametrization space 
analysis to prevent costly computations where they are not needed. For future 
work we propose to attempt to translate the solution from the multi-valued 
to piece-wise linear models, which share the threshold behaviour, but usually 
feature even larger parametrization space and therefore would probably benefit 
even more from reduction techniques. 
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